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Abstract. In this paper, we study the basic problem of counting independent 
sets in a graph and, in particular, the problem of counting antichains in a fi- 
nite poset, from an algebraic perspective. We show that neither independence 
polynomials of bipartite Cohen-Macaulay graphs nor Hilbert series of initial 
ideals of radical zero-dimensional complete intersections ideals, can be evalu- 
ated in polynomial time, unless #P = P. Moreover, we present a family of 
radical zero-dimensional complete intersection ideals Jp associated to a finite 
poset P, for which we describe a universal Grobner basis. This implies that 
the bottleneck in computing the dimension of the quotient by Jp (that is, the 
number of zeros of Jp) using Grobner methods lies in the description of the 
standard monomials. 



1. Introduction 

We approach the basic problem of counting independent sets in a graph and, 
in particular, the problem of counting antichains in a finite partially ordered set, 
from an algebraic perspective. We derive structural considerations and complexity 
results. 

The use of algebraic methods in the study of discrete problems, in particular 
problems in graph theory, was pioneered by Richard Stanley [55] and Laszlo Lovasz 
|20j . from the combinatorics side, and Jiirgen Herzog, Takayuki Hibi, Aron Simis, 
Wolmer Vasconcelos and Rafael Villarreal [17, 25 , 23] from the commutative algebra 
side. The enumeration of independent sets has been approached using Reverse 
Search ([ffi]), the Belief Propagation heuristic ([5]) and Binary Decision Diagrams 
f |27|). to name a few techniques. 

The main algebraic object we will use is the Hilbert Series of the initial mono- 
mial ideal associated with a graph. The problem of computing a Hilbert Series 
is NP-Complete ([I]). There is a standard algorithm (first proposed in [2T]) for 
computing the Hilbert Series of a quotient C[x]/J, where / is a homogeneous ideal 
in C[x]. There are some classes of ideals for which this algorithm finishes in time 
polynomial in the input, e.g. Borel ([T]) and Borel-type ideals ([II])- Open com- 
puter algebra systems (CbGoA [BJ, Singular [T3], Macaulay2 [T2]) implement the 
standard algorithm in subtly different ways. We suggest jT5J Ch. 5] as a general 
reference on Hilbert Series. 

The connection between independent sets and commutative algebra is spear- 
headed by the following construction. 

Definition 1.1. Let G — (V, E) be a graph, with V = {vi, . . . , v n }. The edge ideal 
(H5H22) I' G C C[xi, ...,x n ] of G is defined as 

(1) l' G — (xiXj, for all (vi,Vj) € E). 
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This ideal links independent sets in G and certain monomials. If x a is a monomial 
not in Iq (termed a standard monomial), then it encodes an independent set S of 
G in this way: 

(2) u.eS^^lx". 

This encoding is not one-to-one. For example, the monomials X\ and x\ represent 
the same independent set: {v±}. We introduce a slightly modified version of I' G , 
with which we obtain a bijective encoding. 

Definition 1.2. Let G = (V,E) be a graph. We define the modified edge ideal Iq 
of G as 

(3) I G = I' G + {xl for all Vi ). 

Notice that Iq is zero-dimensional (the origin is the only root), and that its 
standard monomials are square-free, with the degree of a monomial equal to the 
size of the corresponding independent set. The number of independent (or stable) 
sets in G thus coincides with the fc-vector space dimension of the quotient of the 
polynomial ring in n variables over any field k by the ideal Iq. This dimension is 
computed in H2 H2] using the additivity of the Hilbert function in short exact 
sequences. 

In Section [2] we recall the definition of the Hilbert function (see Q) and we ana- 
lyze the instantiation of the standard algorithm for computing the Hilbert Series for 
the ideals Iq- Our main result in this section shows that the recursive calls simply 
correspond to counting independent sets of G that contain a pivot vertex, and those 
that do not contain it. In Section [3j we turn our attention to the subproblem of 
counting the antichains of a finite poset. We present the universal reduced Grobner 
Basis for a family of zero-dimensional radical ideals derived from posets. In Sec- 
tion |4j we specialize our study in the case of Cohen-Macaulay bipartite graphs, 
corresponding to Cohen-Macaulay ideals I' G . Using the characterization in [T7], we 
show that counting independent sets in such graphs is equivalent to evaluating at 2 
the independence polynomial of the comparability graph of a general finite poset. 
Section [5] contains our complexity study. We prove that antichain polynomials 
cannot be evaluated in polynomial time at any non-zero rational number t unless 
P = #P. When combined with the algebraic results from the previous sections 



we deduce Corollaries 5.4 and 5.5 on the intractability of the evaluation of Hilbert 
Series of initial ideals of zero-dimensional complete intersections and independence 
polynomials of Cohen-Macaulay bipartite graphs. We close with a few experimental 
observations in Section |6l 



2. Counting independent sets via the computation of Hilbert Series 

We start by recalling a few definitions concerning Hilbert Series. Let M be 
a positively graded finitely generated C[x]-module (e.g. the quotient C[x]/Iq for 
some graph G). We can write 

M = 0M if 

where M, is the subspace of M of degree i. The Hilbert Function (HF m) of M 
maps i onto dimc(Mi). The Hilbert Series (HSm) of M is the generating function 

(4) HS M {z)=^HF M (i)z i . 

0<i 



If M = C[x]/7 for a monomial ideal I, then HFm{i) is the number of standard 
monomials of degree i (that is, monomials which are not in /). 
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If we take I = Iq for some graph G, as we mentioned in the introduction, 
HF m(*) is then the number of independent sets of size i in G. In this case, the 
Hilbert Series of C[x]//g is a polynomial, called the independence polynomial of G. 
As usual, we denote this polynomial by I(G,x). We refer the reader to [TS] for a 
comprehensive survey of independence polynomias. 

The standard algorithm for computing HS m hinges on the following property. If 
we have a homogeneous exact sequence of finitely generated graded C[x] -modules 

(5) — > M' — > M — > M" — > 0, 
then 

(6) HS M (z) = HS M > (z) + HS M » (z). 

Given a finitely generated graded C[x] -module M and / ^ a homogeneous poly- 
nomial of degree d, we have the following multiplication sequence 

(7) — ► [M/(0 :m (/))](-d) M —> M/fM — ► 0, 

where ip is induced by multiplication by /. Here, (0 :m (/)) ~ {g € M, such that 
gf = 0}, and (— d) induces a degree shift, so that ip is a homogeneous map of degree 
0. Rewriting equation ^ we obtain 

(8) HSm(z) = HS M/fM (z) + z d HS (0:mU)) . 

The polynomial / above is called a pivot. 

Actually, the standard algorithm does not directly compute the Hilbert Series. 
We can see in [HI Theorem 5.2.20] that in the case of the modified edge ideal Iq, 
the Hilbert Series of M = <C[xi, . . . , x^/Iq has the form 

fQl ff9 _HNm(z) 

(9) HSm ~Ji^> 

where HN m{z) is called the Hilbert Numerator. The algorithm computes HN m(z), 
and the series is then obtained by dividing it by (1 — z) n . 

We reproduce the algorithm for computing the Hilbert Numerator of a monomial 
ideal (see [HI Theorem 5.3.7]). 



Algorithm 2.1. Algorithm to compute the Hilbert Numerator of a mono- 
mial ideal I (called HN). 

Require: A set of minimal monomial generators for the ideal /. 

Ensure: The Hilbert Numerator of C[x]/7. 
l: if the minimal generators of / are pairwise coprime then 
2: return J^[| =1 (1 — z di ), where di is the degree of the z-th generator of /. 
3: else 

4: Choose a monomial p as pivot. 
5: fi<-WJ(I:p). 

6: / 2 <-HN(/ + p). 

7: return z d °s(p)/i(z) + / 2 (z). 
8: end if 



Notice that the sets of generators of I' G and Iq described in ([!]) and ^ are min- 
imal. The process of obtaining minimal sets of generators for the ensuing recursive 
calls can be optimized by performing careful interreductions. 

The choice of pivot must satisfy one condition. Namely, 



(10) ^deg(/:p) <^deg(/) and ^deg(I + p) < ^deg(J). 
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Here, ^2 deg(J) denotes the sum of the degrees of all the minimal monomial gen- 
erators of I. Intuitively, this condition says that the recursive calls are made on 
"smaller" ideals, and shows that the algorithm terminates. 

The program GoCoA implements this algorithm, and uses a certain strategy for 
the choice of pivot in step [4] First, it chooses any variable Xi appearing in the most 
number of generators of Iq. Then it picks two random generators containing that 
variable. The pivot is the highest power of Xi that divides both random generators. 

We present a specialized version of Algorithm |2.1| suited for the computation of 
the Hilbert Series of &\x\/Iq for any graph G. 

Theorem 2.2. Let Iq be the modified edge ideal of a graph G. The general al- 
gorithm for computing the Hilbert Series of C[x]/Jg has the specialized version 
presented in Algorithm \2.3\ 

This algorithm has an obvious graphical interpretation. The choice of step\4\ cor- 
responds to choosing a node v of the graph. The recursive calls of step^ correspond 
to counting the independent sets of G that contain v (HScoion) and those that do 
not contain v (HSpi us ). 

Algorithm 2.3. Specialized algorithm to compute the HS of C[x]/7g>. 

Require: The list L of minimal monomial generators of Iq described in Q. 

Ensure: The Hilbert Series of C[x]/J G . 
l: if L consists only of variables and squares of variables then 
2: return (1 + z) k , where k is the number of variables which appear squared 

in L. 
3: else 

4: Choose a variable xi that appears squared in L. 

5: Colon <— a minimal set of monomial generators of ((L) : Xi). 

6: Plus <— a minimal set of monomial generators of (L,Xi). 

7: return zHS Co ion(z) + HS Plus (z) 

8: end if 



Proof. Algorithm |2.3| differs from Algorithm |2.1| in two key steps. In step [T] the 
special version does not check coprimality, as is done in Algorithm |2.1| The other 
difference is in step |4j The specialized version chooses a variable, instead of an 
arbitrary monomial. 

We make a claim that helps us understand why this specialized version is correct. 
In every call to the algorithm, each of the n variables appears in L raised to the 
first or to the second power. Furthermore, in each call, L contains only the powers 
just mentioned and the "edge monomials" XiXj of G such that both Xi and Xj 
appear squared in L. This leads to an obvious graphical interpretation: The list L 
represents the subgraph of G induced by those variables that appear squared in L. 

We prove the correctness of the algorithm by showing that the choice of a pivot 
in Algorithm |2.1| must always yield a variable when applied to a modified edge 
ideal, and that the claim of the previous paragraph is true. 

When the algorithm is originally invoked, every variable appears squared in L. 
Besides the squares of variables, L contains the "edge monomials" XiXj for every 
edge (i,j) of G. This proves that the claim above holds in the first call. 

Assuming that the elements of L have the structure we claim, let us show that 
any choice of pivot yields a variable. Suppose that we employ any conceivable 
strategy for the choice of pivot, always subject to condition (10 1. The pivot p 
cannot be a multiple of any monomial in L. If it is, then Plus = L, and the 
decreasing total degree condition ( 10 ) is not satisfied. The pivot p must then be a 
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product of variables that appear squared in L, but it must not be divisible by any 
"edge monomial." Suppose that the pivot is the product of at least two variables. 
That is, XiXj \ p, where xf and x\ are in L, and Xj,Xj is not in L. Then Plus violates 



the decreasing total degree condition (10 1, because it has the same generators as 



L, plus p. If p = 1, then Colon = L, and this violates the decreasing total degree 
condition. The only valid choice is then p = Xi, for some Xi that appears squared 
in L. 

Once we know that the pivot is always a variable, we can show that the claim 
above holds for Plus and for Colon. In doing so, we also explain the second part of 
the theorem. 

The list of minimal monomial generators for Plus contains all the variables that 
were raised to the first power in L. Furthermore, it must also contain the pivot x^. 
The square of Xi is not in Plus, because Plus is minimal, and the "edge monomials" 
that contained Xi, are not present in Plus. The rest of the generators in L are 
unaffected. Therefore, we have that every variable appears in Plus either squared 
or raised to the first power, as we wanted to show. Plus corresponds to the graph 
obtained by removing the node that corresponds to Xi and all the edges incident 
with it. 

The analysis of Colon is somewhat similar. To obtain a minimal set of monomial 
generators, we just cross out the pivot Xi from every generator in L that contains 
it, and then eliminate multiples. If we had an "edge monomial" XiXj, then Xj 
is in Colon. Therefore, the square of Xj is no longer a generator, and all the 
"edge monomials" containing Xj are also missing from Colon. Again, every variable 
appears either squared or raised to the first power. In this case, we remove the node 
corresponding to Xi, all its adjacent nodes and all the edges incident with Xi or with 
any node adjacent to Dj. 

Let v be the node of G associated with the pivot Xi. The combination step of 
the algorithm reflects the meaning of Colon and Plus: The independent sets of G 
are those of Plus (i.e. those do not that contain v) and those of Colon (i.e. those 
that contain v). 

The algorithm terminates when there are no more "edge monomials" . Since all 
the generators are variables, or squares of variables, then they are pairwise coprime 



and satisfy the stopping criterion of Algorithm 2.1 

A note is in order about the value returned in the base case. Algorithm |2.1| 
returns 

71 

(11) E[(l-**), 

z=l 

where di is the degree of the i-th generator. Since in the specialized case the 



generators are of the form Xi or x\ , expression ( 11 ) has the form 
(12) (l-z)"(l + z) fc , 

where k is the number of variables that appear squared in L. According to for- 
mula ([9]), the value returned by Algorithm 2.3 is the Hilbert Series of C[xi, . . . , x u ]/Ig- 



All these observations show that the graphical interpretation is accurate and 
that the specialized version is indeed correct. 

□ 



3. Partially ordered sets and Grobner Bases 



In this section, we study a family of zero-dimensional radical complete intersec- 
tion polynomial ideals associated with posets, first proposed in 0]. 
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Recall that a poset (or partially ordered set) is a set P, together with a (partial 
order) relation < satisfying 

• a< a, for all a € P. 

• a < b and b < a implies a = b, for all a and b in P. 

• a < b and b < c implies a < c for all a, b and c in P. 

Two elements a and b of P are comparable if a < 6 or if 6 < a. Otherwise, they 
are incomparable. We will usually just write P and drop the partial order relation 
from the notation. 

We can associate to a poset P its comparability graph. 

Definition 3.1. Let P be a poset. The comparability graph G(P) has the set P 
as nodes and there is an edge between two different nodes a and b if and only if a, b 
are comparable in P. 

A subset S of a poset P is an antichain if all the elements of S are pairwise 
incomparable in P. We write A(P) for the set of antichains of P. Note that 
S G A{P) if and only if S is an independent set of G(P). 

Definition 3.2. For any poset P we define the antichain polynomial A{P,x) by 

A(P,x) = I(G(P),x). 

Thus, the k-th coefficient of A{P,x) equals the number of antichains of P with 
k elements and the cardinal |-A(P)| is given by the evaluation A(P, 1). 

Given a finite poset (P, <), we define a polynomial ideal Jp C G[x\, . . . , x n ] by: 

(13) Jp = (xi r, ] [ Xj, for all u< G P). 

Vj<Vi 

Lemma 3.3. Let P be a finite poset. Then the elements of V( Jp) are strings of 
's and 1 's. 

Proof. Let a G V(Jp). Suppose that an element Vi G P is minimal. Then n — x\ G 
Jp, hence is or 1. Now, take any w,, and assume that for every Vj < we 
know that aj is or 1. Note that Xi — Xi Y[ v <v x 3 ~ x i0- ~ Xi Y\v <v x i)- ^ an y 
a,j is 0, then at must be too. If all aj are 1 then ai(l — ai) = 0. □ 

Moreover, we have: 

Theorem 3.4 ( 4 ]). Por any finite poset P, Jp is a radical zero-dimensional ideal. 
Then, it has a finite number of simple zeros. Furthermore, 

(14) \V(Jp)\ = \A(P)\. 

We now show that we can present Jp as a zero-dimensional complete intersection 
by means of generators of lower degree. A standard alternate way of dealing with 
a poset P is to look at the cover relation. Given a and b in P, we say that a -< b 
(read "6 covers a" ) if and only if a < b and there is no c G P such that a < c < b. 
Using this relation we define the ideal 

(15) J P = (xi — Xi Xj, for all u$ G P). 

Lemma 3.5. Let P be a finite poset. Then 

(16) J P = J'p. 

Proof. It is straightforward to see that the varieties of Jp and J' P coincide. We 
show that J'p is radical. Since we already know that Jp is radical, this proves the 
equality. 
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It is enough to prove that the square-free polynomial Xi — xf is in J' P for all 
Vi 6 P. We know this to be true for the minimal elements of P, by the very 
definition of J' p . Suppose we have a non- minimal element Vi in P. Let Vj t , . . . ,Vj r 
be the elements such that Vj k -< Vi, and assume that xj t — an is in Jp for all I. 
First, we observe that x%Xj x — X4 is in J' p for all L Indeed, 



(Xj, - l)(Xi -X? J] 



1 



\ k=l 



Now, consider the following step: 

r 

(xi — Xj J^J ^ife) — ( x i ~ XiXj r )xi J^J Xj k = Xi — x i J j 



1 r-1 
k=l k = l k=l 



Since (xi — xf nl=i ^it) an d (^j — x i x j r ) are m -/pi w e have that Xi — a;f 111=1 ^it 
is also in J' P . If we apply this procedure repeatedly, we eliminate variables from 
the product, and eventually find that Xi — xf is in J' P . □ 



We now take Theorem 3.4 one step further, and give an explicit bijection between 
A(P) and V(J P ). 

Proposition 3.6. Let P be a finite poset. Define the function f : V(Jp) — > A(P) 
by 

f(a) = {v.i 6 P, such that a,i — 1 and a,j — for all Vj > u,}. 
The map f is bijective, and its inverse g : A(P) — > V(Jp) is defined by 

g{S) = a' , where a\ = 1 if 3 vj G S such that v t < Vj and — otherwise. 

Proof. It is clear from the definition of / that no pair of elements of the subset /(a) 
can be comparable for any a € V(Jp), that is, that f(a) is indeed an antichain. 
Reciprocally, let S be an antichain of P and let a' = g(S). We need to see that 
a^(l — a'i Yiv <v a 'k) = f° r an *• This is clear if a\ — 0. When a\ = 1, there exists 
Vj G S with i>j > Vi . By the transitivity of the order relation we deduce that a' k = 1 
for all Vk < Vi and so the equation is satisfied. 

Let a be an element of V(Jp). Let S = f(a) and a' = g(S). We want to show 
that a = a' . Suppose that = 1. Then 3vj £ S such that < Vj, and therefore 
<ii = 1. By a similar argument, if a'i = 0, then ai = 0. □ 

We now describe the universal reduced Grobner basis of Jp — J P . 

Proposition 3.7. The universal, reduced Grobner Basis of Jp is the set Gbp of 
polynomials 

gb i =x\ — Xi Vv, e P, 
gbfj ti ) = XiXj - Xi Vvj < Vi. 



Proof. Lemma 3.3 shows that the elements of V(Jp) are strings of 0's and l's. The 
polynomials xf — Xi are in Gbp, and therefore the elements of V(Gbp) are also 
strings of 0's and l's. Let x = (xi) Vi£ p be a string of 0's and l's. x G V(Jp) if and 
only if V^i G P, (xi = (3 Vj < Vi such that Xj =0)). But this is equivalent to 
x G V(Gbp). Then, Gbp is zero-dimensional, and contains a square-free univariate 
polynomial in each variable (gbi). Therefore it is also radical. This shows that the 
ideal generated by Gbp coincides with Jp. 

We now prove that Gbp is a Grobner Basis for any monomial order <. Recall 
that, given < and a non-zero polynomial p, LT < (p) denotes the largest term of p, 
with respect to <. Clearly, LT < (gb i ) — xf and LT < (gb^ ^) = XiXj. Given any 
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two polynomials in the set Gbp, we show that their ^-polynomial is divisible by 
the polynomials in (Gbp). If we let and q = XkXi — x k be two 

polynomials in Gbp, all possible combinations of the indices i, j, k and £ boil down 
to the following non-trivial possibilities for (p,q) (with i,j,k,£ all different): 

(1) (xf - Xi,x k Xi - x k ) or (xiXj - Xi,x k Xj - x k ) S(p, q) = 0. 

(2) (xj - XuXiXi - xi) =*> S(p, q) = gb, t - gb {Li) . 

(3) (XiXj - Xi,Xix e - Xi) => S(p,q) = gb {hl) - gb {ll) . 

(4) (x^j ~ x i ,x k x l - x k ) S(p,q) = gb {hk) - gb [lk) . 

(5) (x^j - Xi,XjXt - Xj) =>• S(p,q) = gb( t<i) - gb {jti) - 

(6) (xf - Xi, x\ - x k ) or (xj - x k x t - x k ) or (xiXj - x u x k x t - x k ). In all 
three cases, since the leading monomials of p and q are coprime, S(p, q) is 
divisible by (p, q) . 

(7) (xiXj — Xi,XjXi — Xj). This can only hold if i>j < Vj and Vj < Vi, that is, 

Vi = Vj . 

In cases 4 and 5 above, we know that gb ^ m and gb^ e ^, respectively, are in Gbp, 
because a partial order relation is transitive. Therefore, Gbp is a Grobner Basis. 

Finally, none of the polynomials are redundant, all of the leading coefficients 
are one, and the "other" monomial in each polynomial of Gbp has degree 1, so it 
cannot be divisible by any leading monomial of Gbp. Therefore, Gbp is a reduced 
universal Grobner Basis of J p. □ 

We can count the antichains of P by studying J p. We have seen that |^4(P)| = 
|V(Jp)|. It is well-known ([51 Theorem 2.2.10]) that as Jp is radical, it holds that 

\V(J P )\ =dim c (C[x]/Jp). 

The Hilbert Series algorithm could help us to compute dimc(C[x]/ Jp), but it re- 
quires that the ideal Jp be homogeneous, which is not the case. This is circumvented 
by considering an initial ideal of Jp. If < is a monomial order and / is an ideal, 
the initial ideal of LT K (I) is defined by 

LT < (I) = (LT < (p),peI). 

By [7J Chapter 5, Section 3] 

dim c (C[x]/J) = dim c (C[x]/£T<(I)). 

In particular, we have the following equality 

(17) \A(P)\ = dim c (C[x]/LT<(J P )). 

Let Gbp be the (universal) Grobner basis of Jp in the statement of Proposi- 
tion |3.7| By the definition of a Grobner Basis, it holds that for any monomial 
order, LT K (Jp) = (XT< (g),g £ Gbp). Note that this initial ideal has the same 
structure of the ideals Iq in Section [2j In fact, it equals Ig(P)- 

4. Independent sets in bipartite Cohen-Macaulay graphs 

Let G be a graph, and I' G its edge ideal. We say that G is a Cohen-Macaulay 
graph if C[x]//g is a Cohen-Macaulay C[x]-module. The quotient <C[x]/Iq is al- 
ways Cohen-Macaulay, because Iq is zero-dimensional. Cohen-Macaulay rings and 
modules are extensively studied in [TB] , and the article [23] covers Cohen-Macaulay 
graphs. 

Not every graph is Cohen-Macaulay, of course. For example, the path of length 
three (see Figure[T]) has the edge ideal Jp 3 — (xiX2,x 2 x 3 ), defined in C[&i, x 2 , x 3 ]. 
The quotient C[xi, x 2 , £3]/ Jp 3 is not Cohen-Macaulay. It is not even equidimen- 
sional, since the zero set of Jp 3 consists of the plane x 2 = 0, together with the line 
x\ = X3 = 0. 
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Figure 1 . The path of length three P 3 

One particularly interesting subfamily of Cohen-Macaulay graphs are bipartite 
Cohen-Macaulay graphs. 

Definition 4.1. Let G — (Vi U V-2,E) be a bipartite graph. Then G is a Cohen- 
Macaulay graph if and only if G[x]/I' G is a Cohen-Macaulay C[x]-module. 

There is an equivalent characterization, given by the following result. 

Theorem 4.2 ([II]). Let G = (Vi U V 2 ,E) be a bipartite graph. We say that G 
is a Cohen-Macaulay graph if |Vi| — IV2I, and the vertices V\ — {xi, . . . ,x n } and 
V2 = {yi, ...,?/„} can be labeled in such a way that 

(1) (xi,yi) £ E for alli=l,...,n; 

(2) if (xi,yj) £ E, then i < j; 

(3) if(xi,yj) and {xj,yk) are edges, then (xi,yk) is also an edge. 

There are two ways of seeing a bipartite Cohen-Macaulay graph G = (Vi U V2, E) 
as a poset. The obvious way is to set the following partial order on the nodes of 
G: x < y if and only if x = y or x £ Vi, y £ V% and (ar, y) is an edge of G. That is, 
one chooses one of the parts as the "upper" one. 

The other way, which we will consider here, involves a different construction. Let 
G = (Vi UV2,E) be a bipartite Cohen-Macaulay graph. We define a poset Pq as 
follows. The elements of Pq are those of V\. Given x^ and Xj, we set Xi < Xj if and 
only if the edge (xi, yj) is in E. From the transitivity of bipartite Cohen-Macaulay 
graphs, we see that Pq is a poset. 

Conversely, let P be a poset, with elements Xi, . . . ,x r . We build a bipartite 
graph Gp = (V,E) as follows. We set V — Vi U V2, with Vi = {x%, . . . ,x r } and 
V2 = {2/1, ■ ■ ■ ,y r }- We put the edges (Xi,yi) in E for all i, and we have the edge 
(xi,yj) if and only if X{ < Xj in P. In this case, the transitivity of < ensures that 
Gp is a bipartite Cohen-Macaulay graph. 

The following lemma is straightforward. 

Lemma 4.3. The two transformations 

P^Gp and G^f P G 

are inverses. 

We now compare the independence polynomial of a bipartite Cohen-Macaulay 
graph G with the antichain polynomial of the poset Pq. 

Lemma 4.4. Let I{G, x) be the independence polynomial of a bipartite Cohen- 
Macaulay graph G and let A(Pq,x) be the antichain polynomial of its associated 
poset Pq ■ Then 

I(G,x) = A(P G ,2x). 

Proof. The construction outlined above expands every element of the poset Pq into 
a segment in the bipartite Cohen-Macaulay graph G. An antichain S of size k in 
Pg gives rise to 2 k independent sets of size k in the bipartite graph G, since we can 
replace any xt £ S by either the node xt or the node yi of G. It is clear that any 
independent set of G can be seen in this way for a unique antichain S of Pq. □ 
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5. Complexity results 

Is it classically known that it is not possible to count the number of antichains 
of a general poset (that is, to evaluate its antichain polynomial at 1) in polynomial 



time unless P = #P [25]. We extend this result in Theorem 5.3 to the evaluation at 
any non-zero rational number t, by a translation and specialization of [H Theorem 
2.2] to the context of finite posets. We then use our previous results to deduce 
in Corollaries |5.4| and |5.5| the hardness of evaluating the Hilbert function of initial 
ideals of zero-dimensional radical ideals and the independence polynomial of Cohen- 
Macaulay bipartite graphs. 

We start with some definitions. 

Definition 5.1. We define the lexicographic product poset Pi^] of two finite 
posets Pi and Pi as the set Pi x P 2l ordered by the relation (x,i) < (y,j) if 
x < yAx — y =>■ i < j. Similarly, we define the lexicographic product graph Gi[Gy 
of two graphs as the set Gi x G 2 with adjacent to (k,l) iff i is adjacent to k 
or if i — k and j is adjacent to I. 

It is easy to check that -PJPj] is indeed a poset. 

Given a natural number m, denote by K m the poset given by the set {1, . . . , m}, 
ordered with the usual < relation. The associated comparability graph is the com- 
plete graph K m in m nodes, whose independence polynomial equals I(K m ,x) — 
1 + mx. 

It is straightforward to check that the comparability graph of the lexicographic 
product -PJ-P2] of two posets equals the lexicographic product G(Pi)[G(P2)] of the 
respective comparability graphs. We therefore have: 

Lemma 5.2. For any poset P and m G N, the comparability graph of the lexico- 
graphic product P[K m ] equals the lexicographic product of the graphs G(P)[K m ]. 

We are now ready to prove the following theorem: 

Theorem 5.3. Evaluating the antichain polynomial of any finite poset P at any 
non-zero rational number t is -^P-hard. 

Proof. We mimic the arguments in [2] Theorem 2.2]. Suppose, on the contrary, 
that given any poset P on n vertices, there exists an 0(n fc )-algorithm to compute 
A(P, t) for some constant k. Then, given a poset P with n verti ces, consider the 
posets P[K m ] for m = 1, . . . , n + 1. It follows from Lemma 



5.2 



that we can use 

the reasoning in [31 Theorem 2.2] to deduce that that A(P,mt) = A(P[K] m ,t). 
In fact, by [1 Theorem 1]), A(P[K] m ,t) = A(P,A(K m ,t) - 1) = A{P,mt). As 
the posets P[K m ], m = 1, . . . , n + 1 can be constructed in polynomial time from 
the data of P, it would be possible to compute A{P, mt) in polynomial time for 
m = 1, . . . , n + 1. But then, the coefficients ij of A(P, x) — Y]j—q ijX 3 would be 
computed in polynomial time by solving the (n + 1) x [n + 1) linear system with 
invertible matrix M = (My) given by Mjj = i,j = 1, . . . ,n + 1. 

It follows that the number of antichains |-4(-P)| of P would be computable in 
polynomial time by adding y~) . ij . But this counting problem is ^P-complete 
[221. □ 



Combining this complexity results with the algebraic results of the previous 
sections, we have the following two corollaries. 

Corollary 5.4. No algorithm can evaluate the Hilbert Series at a fixed non-zero 
rational number t in polynomial time when applied to initial ideals of radical zero- 
dimensional complete intersections, unless #P = P. 
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Proof. By the results of Sections [2] and |3j the Hilbert Series of the initial ideal 
Ig{p) °f the radical zero-dimensional ideal Jp associated to any poset, equals the 



antichain polynomial A(P,x). The result follows from Theorem 5.3 □ 



Corollary 5.5. There can be no polynomial algorithm to evaluate at any non-zero 
rational number t the independence polynomial of bipartite Cohen- Macaulay graphs 
unless ffP = P. 



Proof. By Lemma 4.4 the independence polynomial I(G,p) of a bipartite Cohen- 
Macaulay graph G and the antichain polynomial A(Pg,x) of its associated poset 
Pg are related by the equality I(G,x) = A(Pg,2x). So, any polynomial algorithm 
to evaluate I(G,t/2) in polynomial time for any Cohen-Macaulay graph G, would 



allow us to evaluate A(P,t) in polynomial time for any poset by Lemma 4.3 The 



result now follows from Theorem 15.31 □ 



6. Some experimental observations 

We tested the three Computer Algebra Systems CbCoA, Singular and Macaulay 
2. The examples we used were the posets consisting of the power set of {1, . . . , n}, 
ordered by inclusion (Boolean lattice). Of the three systems, only CoCoA managed 
to count the antichains for n = 7. These numbers (called Dedekind numbers) 
are known for n up to 13. However, those computations required many hours of 
supercomputer time [T5] , 

The strategy employed by CoCoA for the Hilbert Numerator algorithm seems to 
be generally good. We made some observations about it in [9j. We have also tested 
a recent software package, Edgeldeals ([II]). Edgeldeals allows us to compute the 
Hilbert Series of a modified edge ideal of a graph G by computing the /-vector of the 
simplicial complex associated with the edge ideal of G. The simplicial complex also 
contains a description of the standard monomials of the modified edge ideal of G. 
The computation of both objects (the /-vector and the standard monomials) was 
faster using Edgeldeals for the Boolean lattice, compared to the native Macaulay 
2 implementation of hilbcrtScrics, for n up to 6. 
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